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ABSTRACT 

The leak rates through shaft seals with large prt'ssure drops were simu- 
lated using gaseous hydrogen, or nitrogen flowing through an annulus with a 
nonrotating centerbody. The flows were choked. For concentric or ecienirii 
position of the rotor and parallel or convergent tapered flow passages, data 
and analysis revealed that mass flux or leak rate can be determined 1 rom a 
relation whose normalizing parameters depend on the thermodynamic irilical 
constants of the working fluid and an average flow area expressed in terms 
of the inlet and exit cross-sectional areas. Using thesi- normalized relaiions, 
the flow data for parallel and three convergent tapered shaft seal configura- 
tions are in good agreement. Generalization to any simple gas or gas mixture 
is iniplied and demonstrated in part. 

NOMENCLATURE 

2 

A flow area, cm 

A average flow area, cm 

C normalization constant 

D diameter, cm 

G mass flux, g/cm^-s 

G flow normalizing constant, ‘\/v g/cm‘-s 

G{^ reduced mass flow, G/G* 

h clearance height, cm 

l, length, cm 

m, n exponents 

I’ pressure. Ml’ a 

P|^ reduced pressure, P/P^ 

Rq rotor radius, cm 

T temperature, K 

Tj^ reduced temperature, T/T^. 

compress ibi 1 i ty 


Z 



7 


t density, g/cm^ 

d) mass flow rate, g/s 

Subscripts : 

c thermodynamic critical 

e exit or exhaust 

I inlet 

t throat 

0 stagnation 

I PRODUCT I ON 

A well known rule of thumb in industrial practice is to provide some 
convergence in the design of seal flow passages to minimize the possibility 
of negative stiffness. The "rule" has not been established qualitatively or 
quantatively , however Fleming (1) has found increased stiffness in his analy- 
sis of convergent flow passages with an optimum occurring at an inlet to exit 
clearance ratio of about 1.8. Often overlooked or minimized by large I./D 
are the effects of the inlet. References 2 to 4 demonstrate that in some 
cases, inlet geometry can be of significance in establishing flow rates and 
pressure profiles. Well instrumented 90° sharp edge and Borda type inlet 
results to 105 L/D not only indicate flow separation near the inlet, but 
in some cases of sufficient strength to be classified as critical flow at the 
inlet. For gases the flow recovery to 0.8 Pq was common where Pq is the 
stagnation pressure. 

In many seal designs, the pressure drop across the sealing surface is 
only large enough to choke the flow and some seals are not choked at all; 
however, many current designs require stagnation pressures many times the 
necessary "2:1" for choking. Little data are available on flows througli seals 
with large axial pressure gradients with both concentric and eccentric posi- 
tioning of the rotor in the housing including the effects of convergent coni- 
cal tapers. 

Thus the purpose of these tests will be to determine the mass flux 
through a simulated seal configuration with the fixed "rotor" held in both 
the concentric and fully eccentric positions. 



AFPARAniS AND I NS I'KI 'MKN lAT ION 

rill' bnsii' 1 low facility was of the blowilown type and is desiribed in 
detail in Rel . j. The system was modified somewhat to accomoilati l ho housing', 
which simulated the seal configuration from the space shuttle main luel pump 
interstage seal. A schematic of the modified facility is given as i’ig. I. 
Although the system can handle cryogens as well, the primary objective ol 
these tests is to determine the flow rates for gases nitrogen or liydrogiMi, 
v-/ilh .iiui wi'.hoiii h.!c'<-pri si'.uro . The system bacl<.-;ire:siir-.' was coiUrolUvl bv 
injecting gaseous nitrogen into the exhaust cone. Figure i s a photc'grapli 
of the seal assembly Installed in the facility. A cross-section of the seal 
assembly is shown as Fig. 3. The conical adaptor flanges on each end, while 
not a part of the seal configuration, are necessary to provide propi*r flow 
distribution and measure temperature and pressure. For these tests, the seal 
housing was modified while the centerbody simulating the rotor, remained un- 
altered. A photograph of the instrumented centerbody and the housing prior 
to modifications is shown as Fig. 4. Previous experience indicated that 
minor misalignment could readily be detected in the pressure profi les and in 
the eccentric position, the nonuniform force distribution could be larg.e; so 
in order to minimize misalignment and movement, four tabs at the inlet and 
outlet were used to hold the simulated rotor (centerbody) in a fixed position, 
see Fig. 5. The configurations were machined and data taken in such a way tliat 
two basic housings could provide all the necessary geometric clianges . One was 
machined to give concentric alignment and the other fully eccentric position- 
ing (i.e., contact simulation). Also shown in Fig. 5 are the pressure taps 
whose locations are given in Table 1 along with geometric parameters for the 
simulated rotor and h' ing. It is difficult to see the tapers, yet the (low 
rates are s Ignlf leant .y affected. 

The flow rates were metered using a venturi flowirieter located in the 
bottom of the storage tank and checked against an orifice flowmeter vh.-n ap- 
propriate, that is, when back-pressure control gas was not used. Inl< t stag- 
nation conditions were measured in a mixing ihamber, niit readily visible in 
Fig. 2, md for data reduction the average of fuc pressures within tiu seal 

(‘'02,5,6,d> 



THEOREriCAL BACKGROUNll 

Within a flow passage as shown In Fig. 3, usually one dimensional Mow 


provides an adequate description. Extending some work on two phase choked 
flows and applying the principle of corresponding states (h-lO) one can estah- 
lish choked flow parameters over a large range In temperature, pressure, and 
working fluids. 

The normalizing parameters (b-lO) are G*, Tq, Ff where 


G 


★ 



( 1 ) 


and Fj,, T^,, are the pressure, temperature, density and ci>mpress i b i 1 - 

Ity of the working fluid at Its thermodynamic critical point. See Table 2 
for values of the critical constants for nitrogen, hydrogen, oxygen, argon, 
and methane. Further, it was determined that the choked flow ol real gases 
could be represented over a large rai\ge in pressure and temperature bv 
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where Gg « G/G , Tj^ « 0.33 and 1. < m < 1.02. 

A similar form was advocated by Shapiro fll) tor air. Equilibrium nozzle 
computations using fluids nitrogen, hydrogen, oxygen, argon, and methane 
were carried out using the methods of Refs, b to the most appropriate value 
of the constant appears to be 
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for n = 0.5 and m = 1., with a scatter lor fluid nitrogen of sbqut 3 percent. 
With n “ 0.55 the scatter Is decreased to about I percent with the constant 
Increased to about 0.213. There also appears to be a trend with reduced pres- 
sure such that m ■ 1.02, however tot the extensive charactcrlzat ion ot hy- 
drogen, m - 1.0. Such trends Indicate that for a given rarvge of conditions 
and a given fluid, some improvement over Eq . (3) can be found. 

Recalling that sharp edge Inlets in some cases represent a real flow loss 
of 0.2 Fq, It should follow that experimenlsl flows should be less than 80 
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pfii'ont ot t lu' pti'diiLotl iu>z::lf llows; tiiat ts, in K<( . ( 2 ), tm n 
m =• l.O, tlu- I'onst.'int should bo about O.lb. 

For taporod >;ooinotr ios , tho propor nonn.i 1 i /. i n>; aroa l or mass I lux 
(C = :i'/A) was estimated from somt' experience to be 

2A + A, 

A = Tj ( • ) 

where A^^, is the exhaust flow area and Aj the entrance ) iow area. rh< lal- 
culated average area was found to be 

JA = A^ + A^ t ir(Ry + *’j) * "Rq 

where hj is the clearance height at the inlet and li^, is the clearaiue 
height at the exit, with the rotor radius. Tolerance buildup in both 

the hardware and the data are such that Eq . (S) cannot be ruled out even 
though Eq . (4) will be used herein. 

RESULTS AND DISCUSSION 

In this section we will first discuss the flow data simulation I't normal 
operations (concentric) and then the abnormal- to the point ol rub-case (lullv 
eccentric) for a seal with various convergent tapers. The nominal tapers O.u’, 
0.2^^, 0.3°, 0.4° will be used to characterize the respective geometries. In 
all cases, ambient gases will be used and data taken with back- pressure varied 
from 0.25 to 1.4 MPa. For specifics on the geometries, consult Table I, and 
note that the normalizing area A, given by Eq . (4) is used throughout 

Concentric Geometry 

For the 0.0°, 0.2°, 0.3°, 0.4° taper configurations reduced mass 1 I ux 
data for fluid nitrogen are shown as a function of reduced pressure in Fig. t> 
along with some unpuplished data for the Shuttle seal 1 (unmoilitied version 
ot Fig. 4). In general the results are monotone with a few questionable 
points at high flew rates, which are really at the limit ol the facility. 

Figure 7 represents a similar plot for gaseous hydrogen, except that, 
due to time schedules, no hydrogen data were taken with the 0.2° taper geo;..e- 
try. 

Fully Eccentric Geometry 

The fully eccentric case represents a near-rub and a complex inlet gei'me- 
try subject to separation effects not usually considered in Mow dvnamiis (2-41. 
More will be said about this problem later 
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Figure 8 represents the reduced mass flux results for the 0.0°, 0.2°, 

0.3° taper fully eccentric configurations using fluid nitrogen. Also plotted 
are some unpublished Shuttle seal 1 results and the locus of representing the 
concentric position data. The flow rates are higher than the concentric posi- 
tion by perhaps 6 percent and complex pressure profile behavior was noted. 

The reason may be geometric as strongly suggested from both the nitrogen and 
hydrogen data; however, at this time it remains unresolved. 

Figure 9 is a similar plot for fluid hydrogen with a limited number of 
Shuttle seal 1 points and the locus of representing concentric position data 
is for comparison. No hydrogen data were taken .or the 0.2° taper geometry 
and due to finding of some uncharacteristic and unresolved pressure profile 
results, the 0.4° fully eccentric taper was not machined. 

Normalizing Constant C 

To determine the constant C of Eq. (2), all the data were normalized 
using n >> 1/2 and m > 1, and the resulting value of C was plotted as a 
function of reduced pressure P|^ as Figs. 10 and 11. For the nitrogen and 
hydrogen data there does appear to be a trend with Pq and m ~ 1.02 per 
theory, however the trend is weak and inconclusive. Further, n * 0.5 does 
appear to group these data better than n 0.33, and as noted earlier, 
some Improvement can be made over a limited range with a specific fluid by 
changing the exponents. The nominal values C for these data ere: 

C « 0.153 concentric 

C ■ 0.158 fully eccentric 

This is in reasonable agreement with the theory considering a 0.2 Pg 
pressure loss at the inlet. Under such conditions the constant should be 
C « 0.16. 

Other Factors 

The pressure profiles represent a significant effort, however as noted, 
the results are complex. First the number of pressure taps arc simply not 
sufficient to establish the profiles either axially or circumferentially.. 
Consequently, the pressure profile results are of marginal value but several 
indicative results should be noted for the fully eccentric configurations: 


ORIGINAL PAGE IS 
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1. Thi' overall stiffness of the conical seal C(inf ifinreL i on . as simu- 
lated, appears to be less than or equal to the straiy;ht seal (U.lt" t.ipv r) 
conf igural ion . 

2. It appears that the axial pressure profile tor the 0.0‘’ taper sepa- 
rates within the passage with a subsequent crossover from positive to n«*gativ» 
stiffness. In all cases it appears that the total stiffness is still positiv*-. 
but one must recall that these are static te.sts and such behavior can eite.i 
shaft instabilities under dynamic conditions. 

3. It appears that a "shock" can develop within the seal with the larn'Si 
effect noted in the straight seal. The pressure profiles, the nature ol the 
"shock" and associated separation requires further investigation. 

SUMMARY 

Tlie homogeneous equilibrium choked flow analysis and the principle of 
corresponding states have been extended to the choked flow ol gase.s through 
static configurations which simulate dynamic seal behavior under high pres- 
sure gradients with the following results: 

1. The reduced flow rate for gaseous nitrogen or hydrogen flowing through 
convergent conical or parallel passage geometries can be adequately represented 
by Eq . (2) with m - 1.0 and n = 0.5. It can be inferred using the extended 
theory of corresponding states that any simple gas or mixtures of simple gasat 
will also follow this formula. 

2. For equilibrium nozzle computations, C = 1/5. Therefore assuming a 
0.2 Pq flow loss at the sharp edge the constant (C) would be 0.16. The data 
Indicate (for n = 0.5 and m = 1.0): 

C = 0.153 concentric 

C = 0.158 fully eccentric 

These constants may converge it the minimum stagnation chokini; pressures were 
accurately known for each configuration. 

3. The pressure profiles within these cc’nf igurat ions arc complex and 
separation (shock) effects a change from positive to negative stiffness with- 
in the simulated seal however the overall stiffness appears to be positive. 

4. Assessment of pressure profiles slmul.il ing dynamics requires multi- 
axial and circumferential pressure taps; lurtlu'r work is required here. 
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TABLE 1. - GEOMETRIC PARAMF-ERS 


Concentric 


Diameter 

Diameter 

Diameter 

Nominal 

S 

centerbody, 

Inlet, 

Outlet , 

taper , 


cm 

cm 

cm 

deg 


8.4244 

8.4397 

8.4397 


0.2018 


8.4478 

8.4401 


.2424 


8.4524 



.2625 


8.4587 


mm 

.2861 

Fully eccentric 

8.4244 

8.4409 

8.4409 

0.0 

0.2178 


8.4486 


.2 

.2524 


8.4536 


.3 

.2749 


Seal length = 1.232 cm 


Location of pressure taps 
0.000 cm is at the seal Inlet 


0 ° 

90° 

180° 

270° 

Pq 2 - 2.98 cm 
P05 ■ 0.60 
Pi 0.051 
P 2 0.660 
P 3 0.998 
P 4 1.295 

■ 

Pq^ - 2.98 cm 
Pq 9 - 0.05 
P 5 0.064 

Pg 0.66 

Py 0.988 
Pg 1.295 

P 270 O.bb 
























TABLE 2 


FLOW NORMALIZING PARAMETERS 
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Figure 1. - Schematic of test installation. 
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Figure 2. - Simulated seal assembly installation 




























REDUCED MASS FLUX Gn REDUCED MASS FLUX G 



REDUCED PRESSURE Pr 

Figure 6. - Reduced mass flux as a function of reduced stag- 
nation pressure for gaseous nitrogen flowing through the 
0.0°, 0.2°, O. 3 P. 0.^ taper simulated seal configurations, 
concentric position. 
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REDUCED PRESSURE Pr 

Figure 7. - Reduced mass flux as a function of reduced stag- 
nation pressure for gaseous hydrogen flowing through the 
0.0^. 0. 3^. 0. 4 P taper simulated seal configurations, con- 
centric position. 




REDUCED PRESSURE Pn 

Figure 8. - Reduced mass flux as a function of reduced stag- 
nation pressure for gaseous nitrogen flowing through the 
0.0^, 0.2°. 0. 3° taper simulated seal configurations, fully 
eccentric position. 
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REDUCED PRESSURE Pp 

Figure 9. - Reduced mass flux as a function of reduced stag- 
nation pressure for gaseous hydrogen flowing through the 
O.oP, 0. 3° taper simulated seal configurations, fully 
eccentric position. 
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Figure 10. - Distribution of C with reduced pressure for gas- 
eous nitrogen or hydrogen flowing through simulated seal 
configurations with and without convergent taper, concen- 
tric position. 
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Figure 11. - Distribution of C with reduced pressure for gas- 
eous nitrogen or hydrogen flowing through simulated seal 
configurations with and without convergent taper, fully 
eccentric position. 



